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We perform a parameter study of non-spinning, equal and unequal mass black hole binaries on 
generic, eccentric orbits in numerical relativity. The linear momentum considered ranges from that 
of a circular orbit to ten times that value. We discuss the different manifestations of zoom-whirl 
behavior in the hyperbolic and the elliptic regime. The hyperbolic data set applies to dynamical 
capture scenarios (e.g. in globular clusters). Evolutions in the elliptic regime correspond to possible 
end states of supermassive black hole binaries. We spot zoom-whirl behavior for eccentricities as 
low as e ~ 0.5, i.e. within the expected range of eccentricities in massive black hole binaries from 
galaxy mergers and binaries near galactic centers. The resulting gravitational waveforms reveal a 
rich structure, which will effectively break degeneracies in parameter space improving parameter 
estimation. 

PACS numbers: 04.25.D-, 04.25.dg, 04.70.Bw, 95.30.Sf 



I. INTRODUCTION 

Zoom-whirl orbits arise as a general relativistic phe- 
nomenon of the two-body problem. Such orbits do not 
exist in Newtonian gravity, where the orbits are Kepler's 
conic sections, hence they represent an important facet 
of one of the fundamental problems in general relativity 
(GR). The term zoom- whirl was first used in It 
refers to the orbits of eccentric binaries where tight and 
fast revolutions (the whirls) are separated by phases in 
which the two objects move out to larger distances and 
back in (the zooms). 

The physics behind this effect is precession. For bound 
orbits, one can define the precession per orbit as the angle 
between two consecutive apocenters. In general, preces- 
sion accumulates continuously (with respect to some ex- 
ternal reference frame) and amounts to an excess angle 
beyond the Newtonian motion. Precession is strongest 
for small separations and is therefore significant espe- 
cially for eccentric orbits. In the solar system, the pre- 
cession of the orbit of Mercury due to general relativity 
is 43 arcseconds per century, or (2.9 • 10~ 5 )° per orbit. 
For binary pulsars, precession is not necessarily much 
larger. For the Hulse- Taylor pulsar the precession is 4.2° 
per year, but due to its short orbital period this amounts 
to (3.7 • lfr 3 )° per orbit. Observed are on the order of 
(2 • 10" 2 )° per orbit in some cases [3H5]. A SMBH binary 
model [7] fitted to the optical light curve of the quasar 
OJ-287 predicts (model dependent) orbital parameters 
with precession as large as ~ 40° per orbit. 

In theory, general relativistic orbits with yet larger pre- 
cession can be easily constructed by choosing appropriate 
orbital parameters. This is possible for test particles fol- 
lowing geodesies around a black hole, but also for compa- 
rable mass compact objects in the post-Newtonian (PN) 
approximation [8 . As long as the particle or compact ob- 
ject orbits well outside the innermost stable circular orbit 
(ISCO), the classical picture of a slowly precessing ellipse 



applies. If the orbital parameters are chosen such that 
the object approaches distances close to or even inside the 
ISCO, it may follow an unstable circular orbit for some 
time. After this it either plunges or escapes to larger dis- 
tances (infinity if the motion is unbound), which is the 
zoom-whirl behavior we are interested in. In the whirl 
regime orbits exhibit extreme precession with precession 
angles comparable to or larger than 27r, wrapping the in- 
ner part of the orbit once or even several times around 
its center. 

The basic features of zoom-whirl orbits were first dis- 
cussed in the context of geodesies in a stationary black 
hole spacetime (e.g. [9TfTT]). in extreme mass ratio in- 
spirals (e.g. [DEI El) and PN evolutions QUE]- For 
geodesies, it is a matter of fine-tuning the initial param- 
eters of the orbit to obtain a certain number of whirls. 
In fact, for geodesies the number of whirls can be made 
arbitrarily large since there is no gravitational radiation, 
see |10j for an example with 6 orbits during a whirl. Go- 
ing beyond the test mass limit, including radiation loss 
is a key task, e.g. pQQUQT]. 

The main question about zoom-whirl orbits in full GR 
is how the classic, well known picture of zoom-whirl 
geodesies changes for binaries with comparable masses 
in configurations where radiation damping becomes sig- 
nificant. Naively, we do not expect the binary to radi- 
ate away more than its total mass, i.e. the number of 
orbits is finite since it is limited by the energy and an- 
gular momentum radiated away during each whirl. In 
fact, for comparable masses one might have questioned 
whether it is possible to obtain even a single (full) whirl. 
Since the whirls happen at high velocity and small sep- 
aration (even inside the innermost stable circular orbit), 
the PN approximation is not directly applicable, e.g. [5]. 
However, recently some groups have performed numerical 
evolutions in full general relativity of eccentric black hole 
binaries (BHBs). Zoom- whirl orbits have indeed been 
found, although the number of whirls in these experi- 
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merits is less than three. 

In [TS], Pretorius and Khurana present the first exam- 
ple of a whirl orbit for an equal mass binary. In [191 - 
123] . several examples for the transition from inspiral to 
plunge, radiated energy, angular momentum and the re- 
sulting final spin are investigated. In [23] longer evo- 
lutions of unequal masses and non-vanishing spin with 
up to three elliptic orbits which transition through the 
zoom-whirl regime prior to merger are studied. The no- 
tion of marginally stable circular orbits in background 
spacetimes was shown to be in close resemblance to whirl 
orbits in numerical evolutions of finite mass ratio [T8ll25j. 
The consequences for kicks are addressed in [35]. Impli- 
cations for data analysis are studied in e.g. [2"7H32] . In 
particular [29j [30] point out the potentially deteriorating 
effects in signal processing when eccentricity is ignored 
in the waveform models. 

Eccentric neutron star and mixed binaries in dynam- 
ical spacetime have been studied in |33H36j . and in all 
cases zoom-whirl behavior has been identified. The fo- 
cus in |37H41j is on high-energy collision. Among the key 
results so far is that the total energy radiated can easily 
exceed the 4% of the total mass radiated during the last 
stage of a quasi-circular inspiral. For high-energy colli- 
sion, up to 35 ± 5% have been found [35] • In EH, we 
found at low momentum multiple extrema in the radi- 
ated energy as a function of the initial data, and that 
only a modest amount of fine-tuning is required to spot 
these extrema. These extrema should be compared to 
the variations in the mass and spin of the merger rem- 
nant noted in [23] . 

Choosing different initial data and also different tun- 
ing strategies, these investigations have been performed 
in different regions in parameter space. In the present 
work we focus on an area that has received relatively lit- 
tle attention so far, namely intermediate momenta and 
comparable but not necessarily equal masses. We extend 
the discussion of [35], specifically we consider mass ra- 
tios 1:1, 1:2, and 1:3, and linear momenta that are 1 to 
6, and in one case 10 times the value of a circular orbit, 
although not in all possible combinations. 

Zoom-whirl is sometimes thought to occur beyond a 
certain, rather large eccentricity. Our results instead 
show that whirls can also be found for modest eccen- 
tricities. We give an analysis of the gravitational waves 
(GWs) and how specific features in the radiated energy 
are related to orbital characteristics. 

A prerequisite for zoom-whirl orbits is eccentricity. 
Isolated black hole binaries formed at typical separations 
perform a sufficiently large number of orbits such that 
the orbits become circularized long before entering the 
strong field regime 42J. However, it cannot be expected 
that all binary GW sources are sufficiently isolated and 
hence other effects have to be taken into account. In fact, 
SMBH binaries are expected to be formed in gas/star 
rich environments 43 with potentially large eccentric- 
ities [2U SS]. It is well understood that such binaries 
can gain eccentricity as a consequence of gravitational 



torques exchanged with the circumbinary disk [33] |4"5T - 
49 . Likewise, gravitational interaction with additional 
bodies (Kozai-oscillations, Hill-mechanism, mass segre- 
gation, gravitational focussing, etc.) generically induce 
eccentricity growth on a binary system. Numerous stud- 
ies of such effects [27H29], I50H64] suggest that the eccen- 
tricity of binaries, which emit significant gravitational 
radiation, cannot in general be ignored. Event rate esti- 
mates for eccentric compact object binaries [55] [55] [55J 
suffer from large uncertainties and vary considerably. 
Some studies predict that advanced LIGO should detect 
such sources, but given the large uncertainties this should 
be taken with care. For third generation detectors the de- 
tection range will be larger. Eccentric binary mergers will 
therefore become more interesting sources in the future. 
For supermassive BHBs, pulsar timing arrays will soon 
be able to resolve individual sources in a regime where 
many binaries are still expected to be eccentric (66] . 

Given a population of eccentric binaries, whether 
zoom-whirl orbits are of relevance to gravitational wave 
astronomy depends on several factors. Even if the sig- 
nals are stronger, if excessive fine-tuning is required, then 
the population of strong sources might amount to a very 
small corner of parameter space. Conversely, if little tun- 
ing is involved, then zoom-whirl orbits can be potential 
G W sources even for ground-based detectors [2"7H2"9l 152"] . 
With regard to GWs, for comparable mass binaries with 
astrophysical momenta we loose the unlimited number of 
whirl orbits due to gravitational radiation, but what is 
lost corresponds to a very small part of parameter space 
anyway. In any case, as a matter of principle we should 
be prepared to detect and recognize GWs from all corners 
of parameter space including zoom-whirl orbits. 

The paper is organized as follows. In Sec. [H] we de- 
scribe the basics of our numerical methods, our choice of 
initial configurations, and give error estimates for typical 
runs. We discuss orbital properties in Sec. JHTAI result- 
ing waveforms and radiated energy in |IIIB and phase 



space trajectories in HIC We conclude with Sec. IV 



II. 



NUMERICAL METHODS AND SUMMARY 
OF SIMULATIONS 

A. Method 



We performed a parameter study of the black hole bi- 
nary problem using 3d numerical simulations obtained 
with the BAM code [57-69] . Initial data for black holes 
is computed by the puncture method 70J using a pseudo- 
spectral code [71] . and evolved with the x~ variant of the 
moving-puncture [721 [73] version of the BSSN O [75] 
formulation of the 3+1 Einstein evolution equations. We 
use a 4 th order Runge-Kutta method with 6 th order finite 
differencing. The wave extraction and the calculation 
of the radiated energy is done using a 4 th order accu- 
rate implementation of the Newman-Penrose formalism. 
We extract \?4 and thus also E rac i at extraction radii 



3 



of r GW = 60Af, 80M, 100M, where M the total punc- 
ture mass M (see below). Our grid is a box of typically 
> (640M) 3 size, which is sufficient to keep the boundaries 
causally disconnected from the GWs for most of our runs. 
We employ bitant or quadrant symmetry when possible. 
Usually, the grid consists of 9 levels of mesh refinement 
starting at the coarsest level with resolution of h = 5M 
and increases by factors of two, resulting in the resolu- 
tion of h M/50 at the finest level. The inner, finer 
levels are evolved according to Berger-Oliger timestep- 
ping while the outer levels do not follow the motion of 
the punctures and are evolved at the fixed timestep given 
by the innermost fixed level [SH] US]- The fixed boxes 
have twice as many grid points (for a more accurate wave 
propagation). Since some runs have exceptional settings 
the parameters of our simulations are summarized in Ta- 
ble |lj For our analysis we measured the GW emission and 
radiated energy (normalized to the initial ADM-mass), 
studied the shape of the event horizons, the coordinate 
distance over time and how much time the binary spends 
at a separation D. Moreover we investigated a new way of 
analyzing binary evolutions, namely to look at the phase 
space with the coordinate velocity and the separation of 
the punctures (v,D) serving as generalized coordinates. 
The velocity of the puncture is computed from the shift 
as v = \J f3i(xp)(3 l (x p ), where x p is the coordinate loca- 
tion of the puncture. x p is - as a diagnostic - tracked by 
integrating d t x l p — —f3 l (x p ) using the ICN method as in 
[72]. 

B. Black hole parameters 

The initial data for black hole binaries is characterized 
by a choice of the following parameters. In this work 
we set the spins to zero. Input for the computation of 
the initial data are parameters rrii for the (bare) puncture 
masses, Pi for the momenta, and Xi for the positions. The 
total puncture mass M is defined as M — J2i m i ■ Since 
the global mass scale in vacuum is arbitrary, the masses 
can be characterized by one number, say the symmetric 
mass ratio denoted by v = m\mij (mi + TO2) 2 . 

We choose coordinates in which the punctures are 
initially located on the x-axis, see Fig. [l] For equal 
masses we set x\p = ±D/2 for a coordinate separation 
D. For unequal masses we leave x\ unchanged but set 
X2 = X\m\jm2- For the momenta we choose P\,i = ±P. 
This implies together with the choice of X\^ that initially 
the center of mass is at rest and that mergers happen at 
the origin (except for a small merger kick due to unequal 
masses). Concretely, we consider momenta in the x-y- 
plane given by their magnitude P and an angle such 
that P = (— Pcos0,Psin0,O). Specifying the "shoot- 
ing angle" is equivalent to the choice of an impact 
parameter. The magnitude P of the momenta is chosen 
as a multiple of P qc , which denotes the magnitude of the 
momentum for a quasi-circular inspiral at separation D. 

Given the configuration in Fig. [l] numerical simula- 




FIG. 1. Configuration of initial data. 



tions are parameterized by specific choices for x±, v, P, 
and 0. We set x\ — 10M for all runs, which implies 
D = 20 M for equal masses. For unequal masses, we po- 
sition the larger mass at X2, i.e. m-i > mi and \x^\ < x\. 
Most of the simulations we discuss are for equal masses 
{y = 1/4), but we also consider a few examples for mass 
ratios 1:2 (y = 2/9) and 1:3 (v = 3/16). Following [77] . 
for equal masses at D = 20M the magnitude of the mo- 
mentum for quasi-circular inspiral is P qc = 0.061747A/ . 
We consider P/P qc = 1, 2, . . ., 6, and in one example 
as the most extreme case P/P qc — 10. The direction of 
the momenta is given by € [0,90°]. Here = cor- 
responds to a head-on collision, while for quasi-circular 
inspiraling orbits is slightly smaller than 90° because 
the momentum has a small radial component. The case 
> 90° with initially radially outgoing motion can be 
ignored [75] , 

The ADM mass at the i-th puncture and at infinity is 

M\dm = (1 + u(xi))mi + (1) 

Ma DM = M ADM + M ADM + E bi nd 

= mi + m2 + lim (2ru), (2) 

r— >oo 

respectively, where u is the correction to the conformal 
factor in the puncture framework and Ebind the binding 
energy. Values for M ADM range from 0.994 for lP gc to 
1.2 for 6P gc . Since the momenta are non-zero, we obtain 
larger physical masses M\ DM at the inner asymptotically 
fiat ends of the punctures. The difference between the 
masses m, and M ADAI ranges from 7 • 10 _3 M for lP gc to 
3.5 • 10~ 2 M for 6P qc , and is essentially independent of 0. 

For the main part of this work, we first choose a spe- 
cific mass ratio, in particular we choose between equal 
and unequal masses. Second, we choose one of several 
(low) momentum cases. Third, we vary the shooting 
angle systematically, in particular searching for maxima 
and minima in the total radiated energy, examining the 
number of whirls, etc. There are some obvious alterna- 
tives to set up such parameter scans, say fixing |21j . 
using some measure of eccentricity, the angular momen- 
tum [53], or the binding energy [TH] as parameter. Apart 
from having a simple interpretation as scattering experi- 
ment with fixed momentum size, our setup also describes 
simulations at roughly constant total energy, if in anal- 
ogy to classical point masses the total energy is defined 
as the sum of the kinetic and potential energy (since P 
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TABLE I. List of selected runs with corresponding initial 
data, order of convergence, estimated errors and merger 
times. A 3PN estimate [19, 80 of our lowest eccentric run 
lPqc = 60° gives e ~ 0.5 — 0.6 (depending on which defini- 
tion of e is used) . 



and D are constant while varying 0). Each run amounts 
to 500 - 30000 CPUh (the latter one for lP qc , 6 = 60°), 
which is strongly dependent on how far and how many 
times the orbits zoom out. We implemented Brent's 
method [79] to bracket local extrema in the efficiency 
of converting energy into outgoing gravitational radia- 
tion, for which a small number of runs sufficed. This 
reduced the total number of runs to about a 130 while 
still sampling the parameter space in an adaptive and ac- 
curate way. In retrospect, we found that a golden section 
search 79 is, for the finite accuracy we required, a bet- 
ter choice despite being only first order convergent. The 
parabolic interpolation inside Brent's method chooses the 
new guesses systematically towards the flatter part of the 
asymmetric maxima. 



C. Convergence and error estimates 

We performed a convergence analysis for a represen- 
tative subset of our runs and in general found 4 th order 
convergence in the 22-mode of r^4 and in the radiated 
energy E rac i demonstrating the overall consistency of the 
code with respect to the order of the Runge-Kutta inte- 
grator and the wave extraction routine. The errors due to 
the finite radius of our wave extraction sphere are quanti- 
fied by the deviation from a 1/r fall-off as measured from 
the data taken at three different extraction radii. 

Error estimates based on this analysis are shown in 
Tab. IT and selected convergence plots are shown in 
Figs. [2 and[3j In general, highly eccentric orbits are ac- 
curately treated by the BAM code, and also the presence 
of the rather high momenta considered here can be dealt 
with consistently. The relative errors in the 22-mode and 
the radiated energy both due to finite resolution and ex- 
traction radius is around 1%. For larger momenta the 
error from a finite extraction radius becomes the domi- 
nating error (w 2%). Increasing the initial momentum 
leads to higher amounts of artificial (junk) radiation, en- 
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FIG. 2. Convergence plot of the 22-mode of r*&4 for 
P = lP qc , 9 = 40°. The blue dashed line shows the dif- 
ference of r1le(&4,)\m=2 as computed from the medium and 
low resolution data, and similarly the green (dashed dotted) 
line for high and medium resolution. The solid line rescales 
the latter one assuming 4th order convergence. Both lines 
lie on top of each other demonstrating overall consistency of 
our results. The small noise at the beginning is due to the 
spurious junk radiation of conformally flat initial data. 
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FIG. 3. Convergence plot of the radiated energy for P = 1Pj C , 
O = 40°. The solid red line and the dashed blue one are 
almost on top of each other, demonstrating a convergence 
order very close to 4 consistent with our numerical scheme. 



hances the ADM mass of the initial time slice and reduces 
the BH horizons. At some point these effects contami- 
nate the solution in the sense that its physical relevance 
becomes questionable. However, we limited our data set 
to those regimes where the artificial radiation is either 
entirely negligible or at least small in comparison to the 
physical radiation. The P = 10P gc sequence represents 
an exception, but we only used it in the context of the 
Hawking limit and as an approximate extrapolation of 
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our data set to the ones obtained by other groups at 
larger momenta. 

The error due to the junk radiation arising from the 
conformally flat initial data can be reduced to some ex- 
tent by choosing a sufficiently large initial separation. In 
the very high momentum case we needed both large sepa- 
rations and also much larger resolution until the radiated 
energy results converged, but once the appropriate res- 
olution is used the accuracy compares favourably with 
the other results. For yet higher momentum runs we re- 
fer to [37J , who study momenta beyond the rest-mass 
dominated regime. 

Convergence for unequal masses can be shown only at 
the higher resolutions. For the resolution we used the 
1 = 2, m = 2 mode converges at second order - a com- 
mon tendency when being at the edge of the convergent 
regime. The somewhat lower accuracy for larger mass ra- 
tios is a well-known effect of the gamma-driver condition 
we use (rj = const) . In [5TH55] it is shown that a general- 
ization of this condition (77 dependent on the local mass) 
leads to an improvement in accuracy. 

Analyzing the dependence on resolution shows that the 
derived errors in the energy are not behaving according 
to a Gaussian distribution. There is a skewness in the 
actual (unknown) distribution of our measurements such 
that higher resolutions systematically produce higher en- 
ergies. Hence, our errorbars should be slightly more ex- 
tended towards larger values of E rac i- 

Summarizing, the simulations presented here do not 
pose new challenges to the numerical scheme, although 
there are specific requirements for accuracy in the pres- 
ence of whirls together with long runtimes. In these cases 
there is a high sensitivity to the parameters and during 
the long evolutions numerical errors accumulate. Never- 
theless, those evolutions have similar convergence behav- 
ior and error estimates, and only require a higher resolu- 
tion to obtain convergence. 



III. RESULTS 
A. Orbital properties 
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FIG. 4. P = P qc , 9 = 42°. From 9 = 0° to about 40°, 
the black holes collide within 200A/ of evolution time and 
perform less than one orbit before merger. The waves show 
a merger signal with a brief ring-up and characteristic ring- 
down. In the limit of the head-on collision, 0—^0°, the 
22-mode vanishes and the 20-mode becomes the dominant 
mode. 
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FIG. 5. P = P qc , 9 = 48°. There is about one full whirl 
in a range of ±1° around this shooting angle followed by the 
merger. The waveform clearly shows a wave associated with 
the whirl. Its amplitude is smaller than that of the ensuing 
merger signal. The diameter of the whirl is smaller than the 
innermost stable circular orbit of a Schwarzschild BH with 
the same total ADM mass. 



To prime the discussion of the orbits, we first con- 
sider several examples of puncture tracks for equal mass 
binaries with P = P qc and P = 5P qc , see Figs. [4] - 



10 It is helpful to read the captions of these figures 
in sequence. Shown are the puncture tracks in the x-y- 
coordinate plane in the upper panels, and the 22-mode 
of the waveforms in the lower panels. The waveforms are 



further discussed in Sec. IIIB The figures show two se- 
quences of runs for two momenta that explore how the or- 
bits change when the shooting angle is varied from small 
to large. 



1. Classification of orbits 

For any choice of mass ratio v and initial separation 
D, we can in principle fill in a "phase-diagram" as shown 
in Fig. [II] which labels orbits in a P-0 plot. The main 
classification is whether initial parameters P and O lead 
to orbits that are bound (implying capture and merger) 
or unbound (escape to infinity). In Newtonian gravity, 
we only have to check whether the kinetic energy exceeds 
the potential energy, or equivalently whether the binding 
energy is positive or negative. In general relativity, this 
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FIG. 6. P — Pqc, O = 50°. For a shooting angle two degrees 
larger than that leading to the strong whirl, there is a close 
encounter with a precession of about half an orbit, followed 
by a zoom out to about three times the radius at pericenter, 
followed by a short inspiral and merger that starts with signifi- 
cantly reduced eccentricity. Note the comparatively small and 
short wave pulse associated with the close encounter, again 
at about 200Af of evolution. 



FIG. 8. P = 5P qcy 6 = 14.15°. Example for initial momen- 
tum which is significantly larger than that of quasi-circular 
orbits, and which can easily produce unbound orbits. Zoom- 
whirl orbits are found for much smaller shooting angles than 
in Fig. [4] -[7] There is one whirl, and a short zoom followed 
by a merger. Due to the additional kinetic energy, the whirl 
signal increases in amplitude and exceeds the merger signal. 
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FIG. 7. P = Pqc, <d = 60°. Increasing the shooting an- 
gle beyond 50°, one can find an increasing number of elliptic 
orbits. Early on the orbit resembles the classical picture of 
a (strongly) precessing ellipse. The plot shows a transition 
through plunge through a full whirl phase at the onset of 
merger with a clear corresponding wave signal. 



FIG. 9. P = 5P qc , 8 = 14.20°. The larger the momen- 
tum, the more sensitive the orbit becomes to the choice of 
the shooting angle. A small change in angle compared to 
Fig. [8] leads to a much larger zoom out to an apocenter dis- 
tance of 12M before it merges at the next encounter. The 
initial whirl, however, is almost unchanged, highlighting the 
analogy to unstable circular orbits. 



distinction is sometimes only possible a posteriori since 
the gravitational waves and the associated loss of energy 
and angular momentum are only known after the Ein- 
stein equations have been solved. Solutions to the evolu- 
tion problem define the dividing line P = Pb u {v, D : Q) in 
Fig. [11] Orbits with P > P bu are unbound, orbits with 
P < Pbu are bound. 

A simplified, a priori upper limit on the momentum P 
that ensures boundedness is Pb := Pb u {v,D,0 = 180°), 
which is independent of 0, cmp. Fig. If the mo- 
mentum P does not suffice to escape in the direction 



= 180° (for which radiation losses are minimized), 
then the orbits are bound for all 0. Here we use the 
assumption that the black holes are not spinning. Ap- 
proximating the minimal radiation loss in the "head-off" 
direction by zero, we compute a simple estimate of Pb 
based on the binding energy in Fixing 6 = 180° 
and D = 20M we iteratively compute initial data with 
varying P to obtain the binding energy Ebind- P is then 
defined as P := P(E bind = 0,6 = 180°) resulting in 
P = 0.085(4 ± 3)M. There is some error since we end 
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FIG. 10. 



5Pqc, 



6 = 14.30°. Enlarging the shooting 
angle further compared to Fig.[9]results in a full whirl followed 
by a zoom to infinity (unbound orbit, no merger). 
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FIG. 11. This plot sketches the end state of eccentric black 
hole binaries in the plane spanned by our parameter choice 
for the initial data. The evolutions located in the grey shaded 
regions can be judged bound solely based on the initial data. 



the iteration at some point and since radiation effects 
are ignored. For example, the momentum of the quasi- 
circular orbit leads to bound orbits for all angles since 
P qc <P^ 1.377/V 

For brevity, we will refer to the set of configurations 
satisfying P < P as the elliptic regime. On the other 
hand, orbits with P > P form the hyperbolic class. Note 
that this terminology skips over the fact that orbits in the 
hyperbolic regime may still lead to a merger provided O 
is small enough. 

Inverting Pb u (u,D,Q), we define @bu(v,D,P), the 
shooting angle between bound and unbound orbits, as 
a function of P. An a priori upper limit for bound orbits 



is given by 6 ffeom = Q bu (v,D,P = oo), see Fig. [II] In 
practice it is tricky to study large P due to limitations in 
the construction of initial data. Conceptually, however, 
we can think of this limit as a geometric constraint based 
on the finite size of the black holes, i.e. the idea is that 
the two black holes must merge when their event hori- 
zons touch. Using Euclidean geometry, Q geom is given by 



sin(6 9eoro ) 



r/D, where d n 



is the separation 



of the punctures at the time of the merger. However, 
the size of the black holes depends on the gauge. The 
Schwarzschild radius for a mass m is 2m in Schwarzschild 
coordinates, m/2 for isotropic coordinates, and depend- 
ing on the moving puncture gauge somewhere in between 
for the numerical evolutions. We therefore use the numer- 
ical result for d. 



merger • 



For equal masses we find that a 
common event horizon appears at a coordinate distance 
of about d merger « 1.76 to 1.95M (with a slight drift 
towards smaller values with increasing momentum). For 
an initial separation of D = 20M This estimate leads to 
a geometric limit of @ geom = 10.5° using the Euclidean 
formula. This limit does not appear to be very restric- 
tive for low momenta, but it is not in contradiction to 
the runs of this study, either. All our simulations with 
< Qgeom end in a merger. 

The determination of the ultimate fate of a system out- 
side the above ranges requires a full numerical evolution. 
Here a bound system can be defined by the (future) for- 
mation of a single event horizon, which is expensive to 
compute numerically. In our evolutions we use a crite- 
rion on the lapse at the center of our grid to determine a 
merger time. We justify this approach by a direct com- 
parison with an event horizon finder 84, 85] . The merger 
time t m is approximated by t m , the time by which the 
lapse at the center of our grid has dropped below a = 0.3. 
This is near the analytical value of a single Schwarzschild 
black hole in the same and similar gauges [8BTI88) . We 
have chosen a moderately long evolution among the ellip- 
tic category and get t m = 484.175M and C = 485.524, 
accurate to within At m /t m — 0.0028. We should men- 
tion, however, that the lapse criterion gives worse an- 
swers when the punctures move too fast, because the 
value a = 0.3 is motivated by a Schwarzschild spacetime 
and hence is not well adapted to a boosted black hole. 
We used the lapse criterion to estimate the merger times 
and list them in Tab. HI Those values are also used in 
Figs. [19] and [20] 

Even if one performs a numerical evolution it can be 
difficult to determine whether an orbit is unbound. The 
absence of a common horizon is only a necessary but 
not a sufficient condition for unboundedness. If a merger 
does not occur after a given finite time, the question is 
for how long the simulation has to be continued to settle 
whether the binary is bound or unbound, and in principle 
this time can be infinite. A practical, approximate crite- 
rion can be given in terms of the initial binding energy 
Ebind and the energy radiated in GWs (see Sec. IIIB) 
during the first encounter. Without gravitational radi- 
ation Ebi n( i is a constant of motion and the orbits are 



find, unsurprisingly, that all orbits with Eund < also 
merge in our evolutions. We judge an orbit to be dynam- 
ically captured when the energy radiated during the first 
encounter exceeds the initial (positive) binding energy. 
This shortens the runtime to determine whether a run is 
unbound significantly because we do not have to track 
the black holes to larger and larger distances. Such a cri- 
terion is applicable close to the threshold between bound 
and unbound runs, although a few marginally bound runs 
may be incorrectly labeled unbound (but runs are labeled 
bound correctly). 

We conclude with remarks on the relation to peri- 
odic orbits. Within the category of bound orbits there 
is a detailed classification scheme based on periodic or- 
bits which is complete when neglecting radiation effects. 
In this classification [UJ EH] one indexes all closed or- 
bits with a triplet of integers (z,w,v), where z is the 
number of zooms within an approximate 2ir period (i.e. 
the number of "leaves"), v is the stride over the leaves 
(1 < v < z— 1), and w is the number of whirls. The total 
precession angle is 2ir(w + -). The question is whether 
this classification still works in an approximate sense for 
BHBs with radiation effects. Especially near the merger 
of comparable mass BHBs, the orbits shrink significantly 
and may not be well represented by a single periodic 
orbit, but rather by a sequence of them. Our findings 
imply that the longest whirls associated with the largest 
precession angles (largest w) occur for momenta with P 
slightly larger than P and are very close to a precession of 
2ir. We also find that the dependence on P is weak and 
beyond P > 2P qc compatible with the statement that 
it only depends on the mass ratio. Radiation damping 
seems to limit the length of the whirl phase for larger P, 
although there may be artifacts due to the initial data. In 
terms of periodic tables this means that we typically find 
preferred subsets of periodic orbits that approximate our 
evolutions best. The number of whirl-orbits w is clearly 
limited by the efficiency of gravitational radiation. For 
equal masses w — 1 seems to be the largest w one can 
obtain. For larger mass ratios w — 2 should also be- 
come possible somewhere beyond a mass ratio of 1:3. In 
the regime we are probing orbits with z = 2, z = 3 and 
v = 1 are favored. However, our data set contains too 
few data points on different mass ratios to make a strong 
statement. 



2. Examples for orbital dynamics of BHBs 

We describe the main aspects of the orbital dynamics 
that we find in our data set using the categorization in- 
troduced in the previous section. First we consider equal 
mass BHBs in the elliptic regime. All equal mass runs 
start at D = 20M (P = 10P qc has D = 50M) in such 
a way that D shrinks. Obviously, the ensuing evolution 
depends on the values of P and O. 

We discuss the orbital dynamics from low to high O 
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FIG. 12. Coordinate distance in the P = lP qc sequence as 
a function of time. For larger O more eccentric orbits in the 
inspiral become possible before merger. We simulated up to 
6(+l) orbits for O = 60°. Note the plateaus due to whirls 
before merger in the 48°, 51°, 60° cases in contrast to the 
short plunges in the other cases. 



for P 



P qc . Puncture tracks for some values of O are 



shown in Fig. [4]-[7j while Fig. fl2] shows the coordinate 
distance D. The insets of Fig. |14| show puncture tracks 
for some additional values of 6. 

At low O (or equivalently for high eccentricities) D 
monotically shrinks leading to a rather prompt merger 
without completing a single orbit. The runs with larger 
have correspondingly higher initial orbital angular mo- 
mentum and manage to resist the strong gravitational 
pull for longer so that the merger time steadily grows. 
For O w 46° the punctures complete one orbit before 
merger. At yet larger O > 48° the orbits begin to ex- 
hibit a circular phase (the whirl) which is maintained for 
longer as O is increased. However, at O ss 48.5° the orbit 
leaves the circle again towards larger radii (the zoom) de- 
laying the merger significantly. In this range of O there 
is high sensitivity to the initial data (concerning merger 
time as a function of O). A mild increase in O leads to a 
much larger t m because the BHs slow down as they move 
out before falling back. In the limit O — » 90° the peri- 
center passages become shorter while the apocenters and 
pericenters become increasingly degenerate. The peri- 
center moves out with O, hence the BHs do not cross 
their mutual gravitational potential as deeply and conse- 
quently not as much radiation occurs, enabling more and 
more orbits before merger. 

Concerning the amount of precession, we see that al- 
though our evolutions start somewhere beyond the apoc- 
enter (e.g. Fig. 14), the orbits exhibit a huge precession 
of roughly tt and close to 2n for O = 48.5 (followed by a 
tiny zoom). Even for the smallest eccentricity we stud- 
ied O = 60° (e ~ 0.5) one can see that the ellipses still 
have precessions as large as 27r/3, meaning that over the 
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course of the whole evolution the accumulated precession 
amounts to more than two entire orbits. Those values by 
far exceed the amounts of precession known from mildly 
relativistic systems like the famous Hulse- Taylor pulsar 
[3] with a precession of q = 0.0037° per orbit or the bi- 
nary pulsar [4] with q = 0.0044° per orbit. Furthermore 
for these systems and their correspondingly milder grav- 
itational radiation there is precession not only from one 
apocenter to the next but also precession of the multi-leaf 
clover as a whole very similar to the findings in studies of 
periodic orbits QTJQ3] (or nearby aperiodic orbits). This 
additional, peculiar precession effect is indeed small for 
our evolutions as well, though not negligible. 

Switching to the hyperbolic class the additional pos- 
sibility arises that the BHs just fly past each other, de- 
flecting their trajectories and escaping to infinity. This 
gives rise to the merger / fly-by threshold (see insets in 
Fig. 15 ) which we will discuss later. 



We again describe the orbital phenomenology from 
low to high values of 0. The qualitative features of 
low evolutions are the same as in the elliptic cate- 
gory. The actual values of that lead to analogous fea- 
tures/characteristics (one complete orbit, a whirl, maxi- 
mum in E ra( i, etc. ) decreases for increasing P. This is 
in agreement with the expectation that for larger P one 
has to shoot the BHs closer to each other compared to 
the corresponding lower P evolutions in order to obtain 
a qualitatively similar behavior. Again as in the elliptic 
category for larger values of the whirl phase is fol- 
lowed by a zoom. Depending on P there now is a finite 
range in where the BHs do not escape to infinity, but 
reach an outer turning point (like for elliptic orbits) and 
fall back ending in a delayed merger. Beyond a certain 
(momentum dependent) value of M {, the BHs are simply 
deflected or fly-by each other. 

Inside the hyperbolic category, in our simulations we 
only find orbits exhibiting whirls during the first en- 
counter and never thereafter. This behavior agrees with 
the interpretation that too much angular momentum is 
radiated during the first whirl to have another whirl 
episode. Another way of explaining this is to realize 
that the first pericenter distance during which dynam- 
ical capture occurs is already within any (quasi-) stable 
orbit. On the next encounter the binary will have lost 
additional angular momentum and will have a yet smaller 
pericenter separation. The system therefore is likely to 
merge on the next encounter. It is unclear whether for 
high momenta there can also be cases where after a first 
whirl the orbital parameters fall into the narrow window 
for a second whirl. 

In geodesic motion there exist solutions that escape to 
infinity after a full 2n whirl. Like in previous studies, 
we were not able to find such orbits in the comparable 
mass case. This is most probably due to excessive loss of 
energy and angular momentum during the whirl. It is an 
open question whether this statement holds for general 
momentum P, but if it occurs then for momenta lP qc < 
P < 2P qc or P > 6P qc . 



We proceed by analyzing precession effects and discuss 
resemblances to periodic orbits. For a given P the pre- 
cession angle shrinks with increasing when approach- 
ing the threshold as expected. The maximal amount 
of precession we find is slightly larger than in the ellip- 
tic category. We clearly recognize patterns known from 
periodic orbits. For the P — 2P qc sequence we find 
z = 2, z = 3, z = 4 orbits. The main difference to pe- 
riodic orbits is that the orbits end in a merger after the 
first leave has been traversed because of the severe radi- 
ation losses. For instance P — 2P qc , = 25.1° resembles 
the z — 3,v — l,w = orbit with q = 2n/3. When 
decreasing by small amounts, the resulting orbits typ- 
ically show the same amount of precession (only D per 
shrinks with 0). At some point there is a transition to 
another multi-leaf clover and the precession amounts to 
a value of q = it and is now similar to the periodic orbit 
labeled z = 2, w = 0,v = 1. 



3. Unequal mass BHB and geodesic limit 

Next we extend the discussion to unequal mass BHBs. 
By doing so we move towards a region in parameter space 
which can be increasingly well described by geodesies. In 
fact, it has been in the latter regime where zoom-whirl 
behavior was studied first [T]. This begs the following 
question: Given a binary at a finite mass ratio, how far 
away is it from the geodesic limit? 

The fact that zoom whirls can be found not only for 
geodesies, but also for equal masses suggests that zoom 
whirls also occur for intermediate mass ratios and adds 
to their expected astrophysical relevance. Indeed we 
can confirm (see also [Ml 131]) the presence of zoom- 
whirl behavior for mass ratios 1:2 (v w 0.2222) and 1:3 
(y = 0.1875) (see Fig. [24]. 



As the mass ratio departs from unity, gravitational 
radiation decreases, e.g. [5S], which is consistent with 
the trend to the geodesic limit. In the eccentric case 
we find that qualitatively a similar statement still holds. 
We point out, though, that there is a non-trivial depen- 
dence on (or inverse eccentricity). In particular, the 
maximum in E ra d/MADM (see Sec. IIIB2 and Fig. 18) 



is close to the equal mass values for the mass ratios we 
have probed. 

For lower symmetric mass ratio v we do not find sig- 
nificantly longer whirl phases in our data sets. It is to be 
expected of course that for some mass ratio beyond 1:3 
the whirl phases eventually will be longer and asymptote 
to the geodesic limit. Highly eccentric binaries with mass 
ratios up to 1:3 are in this sense still far away from the 
geodesic limit. 

We find evidence for the analogy of zoom-whirl dy- 
namics and unstable circular orbits by investigating the 
orbital radius during the whirl phase for various configu- 
rations. Consistently, the whirl radius decreases with in- 
creasing P, which we will refer to as the tightening of the 
whirl. This is consistent with earlier studies [T51 |2"TI |2"3"] . 
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in which it was found that the spin of the merger remnant 
increases with the initial angular momentum parameter, 
which implies a smaller radius for the unstable circular 
orbits. 

Next, we investigate geodesies to derive lower limits 
on the shooting angle that separates merging from non- 
merging evolutions, Qbu- (The corresponding values from 
our evolutions can be seen as a vertical dividing line in 
Figs. [15) [17J and [18]) Since we find that 0{, u decreases 
monotonically with increasing P, we expect this lower 
limit to be most restrictive for large P. The idea is 
analogous to the capture/escape cavities for a photon in 
Schwarzschild spacetime in [50] . 

A null geodesic in the Schwarzschild spacetime on a 
circular orbit is located at a radius equal to the so-called 
photon orbit r photon = 3m, which leads to the limit: 
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The proper computation of a null geodesic in 
Schwarzschild spacetime [90] leads to 
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FIG. 13. Higher modes of rlZe^i) summed over m for the 
equal mass case. Clearly one can see that I = 2 modes are 
- just as in the quasi-circular case - the dominant contribu- 
tion. The plot also reveals that higher modes exhibit a much 
more significant contribution from higher Z-modes. All odd I 
vanish within numerical error as expected from the quadrant 
symmetry of equal mass, non-spinning BHB. 



The same calculation for marginally bound circular orbits 
yields 21.8°. Fig. 17 indicates that neither of these lim- 



its apply to our evolutions, because there are unbound 



orbits with < 9 



geod 
b u 



Clearly, the assumption of a 



Schwarzschild spacetime is not a good one. 

From [211 I23j we know that the merger remnant in 
our settings will settle down to a Kerr solution with 
spin parameters between 0.6 < a < 0.823 with only 
weak dependence on the initial conditions. Despite the 
fact that the Kerr metric docs not describe the space- 
time at merger, it may be a better approximation than 
Schwarzschild. The same estimate as above for Kerr 



0.6, QIT 



11.4 



spacetime yields Q? = 14.9° for a 
for a = 0.0823 and Q 9 b e u ° d = 5.7° for a = 1. The 6-values 
for 0.6 < a < 0.823 correspond rather well to the shoot- 
ing angles separating bound from unbound runs in the 
higher momentum cases despite the fact that the limit 
from null-geodesics to (finite-size) equal mass binaries 
is by no means straightforward. We will use this anal- 
ogy in interpreting our results on the radiated energy in 
Sec. |IIIB 2| based on the tightening of the whirl orbits 
associated with a larger spin of the merger remnant. 



B. Radiation properties 

1. Waveforms 

The methods used to compute quantities characteriz- 
ing the GW content of the spacetime are described in 
e.g. |67j . Here we demonstrate how the orbital dynamics 
as described in Sec. |III| 4. are reflected in the GW signals. 



The waveforms of quasi-circular binaries are rather well 
understood. To a certain extent merger waveforms as 
they arise from evolving quasi-circular binaries can be 
very similar to the ones seen in low eccentricity evolutions 
provided the binary circularizes before merger. For large 
eccentricities it is, however, natural to expect deviations 
from a quasi-circular BHB. We observe differences in the 
waveforms throughout the evolution including inspiral, 
onset of merger, coalescence and ring-down. Any im- 
prints left from the eccentric inspiral have to be radiated 
away during this process, because the final spacetime can 
be described by the Kerr metric. In fact, the merger 
remnant reveals a different signal during ring-down |22) . 
In particular, quite generically high eccentricity is corre- 
lated with an amplified ring-down signal. 

The inspiral features show some level of agreement 
with |91j and PN models for such waveforms are known 
analytically to 2PN order [S2] (see the first compari- 
son between numerical waveforms and Post-Newtonian 
ones in the eccentric regime [93]). However, features as- 
sociated with zoom-whirl behavior (see Fig. [7] prior to 
merger) are exclusive to the strong field and thus have 
to be dealt with using the tools of numerical relativity. 
These inspiral signals will be observable by future GW in- 
terferometers such as LISA [SH[S5J or eLISA/NGO [SB], 
DECIGO [ST] or the ET-telescope [SHUSSJ. 

We discuss typical waveforms of a representative subset 
of our evolutions. It is illustrative to go through Figs. [4]- 
T0] 14 and their captions. Our main focus is on the rich- 
ness in information stored in eccentric BHB waveforms in 
contrast to quasi-circular ones because of the promising 
implications for data analysis, see [32) . 
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Already the 22-mode shows obvious differences which 
become larger in other modes. For example the I = 2 
rn = 0-mode of a quasi-circular orbit looks just like a 
smaller amplitude version of the I — 2 m — 2-mode. 
In the eccentric case they contain completely different 
features. We plot the higher I modes summed over m in 
Fig. [13] for the equal mass case, for which (without BH 
spin) only even Z-modes contribute by symmetry. We 
have computed the I < 8 modes and find that the I = 2 
is still the largest contribution, but I — 4 has a significant 
contribution throughout the merger and I = 6 and I = 8 
close to the maximum. 

As an example of how different waveforms of binaries 
with eccentricity and mass ratio away from unity can be, 



we show in Fig. 24 the waveform and orbital trajecto- 
ries for the mass ratio 1:3. Clearly the features induced 
by eccentric unequal mass BHB give rise to waveforms 
which effectively break degeneracies in parameter space 
[321 [TOO]. 

In this work we do not construct waveform templates. 
Longer runs will be needed in order to achieve a match 
to a PN waveform because of the small separations at 
pericenter. Performing wave extraction at larger radii is 
also clearly desirable in this context. With current codes 
this could be done at an acceptable computational cost. 



2. Radiated energy 

We compute the energy E rad radiated away in GWs 
and analyze these results together with the orbital dy- 
namics. For the elliptic orbits we add an estimate of the 
radiated energy of the past evolution ~ ~Eu n d{t = 

0) to E ra( i. Using this estimate we implicitly assume that 
the binary was isolated in its entire past. The actual 
value Eu n( i(t — 0) for the P = P qc sequence turns out to 

be Kad ~ ~ E bind(t = 0) fa 0.0057 ± 0.0001. We nor- 
malize E ra d by the ADM-mass of the initial time slice, 
Madm (t = 0). The resulting quantity is what we call 
the "efficiency" of gravitational radiation. 

The results of all our evolutions are presented in Figs. 

The different lines (colors, sym- 



15 16 



17 and 24 



14[ 

bols) in these plots correspond to different initial mo- 
menta and each line shows the efficiency of gravitational 
radiation as a function of 0. 

The first global feature to notice is that gravitational 
radiation becomes much more efficient for higher mo- 
menta. We give the maximal efficiency for 7 initial mo- 
menta P. So far the largest value 35 ± 5% was reported 
in pi|J EH [3jH 1102] (in which the punctures have coor- 
dinate velocities of v = 0.94). In our data set we come 
rather close to this limit (see Fig. 16 1. The challenge in 



these studies arises from the growing significance of un- 
physical radiation content that is associated with the con- 
struction of initial data. Here we did not intend to push 
this limit further, but this shows that we have probed 
part of the parameter space close to the limits of former 



investigations. As we shall demonstrate (see Fig. 17) the 



sampling is quite exhaustive and allows us to probe zoom- 
whirl behavior in a large class of orbits. In particular, one 
of our important findings is model P = 1P QC , = 60° (el- 
liptic class) with several close encounters before merger, 
see Fig. [7J The initial eccentricity is as low as e ~ 0.5. 
This is a value within typical estimates of supermassive 
BHBs in galaxy merger scenarios following star- or disk- 
driven hardening |49| and also a value found for inspiral- 
ing binaries near galactic cores [53] which are driven to 
very similar eccentricities via the Kozai mechanism. 

For low momenta and the mass ratios under consider- 
ation the shooting angles for the largest number of or- 
bits in general neither coincide with the maxima in E ra d 
nor do they coincide with the unstable, circular (whirl- 
like) orbits merging right thereafter. Generally the maxi- 
mum in E ra( i inside the hyperbolic regime lies close to the 
merger/fly-by threshold. However, in the limit P — » 00 
there appears to be a growing amount of degeneracy: the 
unstable-circular orbits actually seem to coincide with 
the most efficient radiators. In the next section we will 
give an interpretation for this behavior. 

For low we find, in agreement with previous studies 
[15] , that the radiated energy quickly drops to the small 
amounts known from head-on collisions [30] - This drop 
can clearly be seen for every initial momentum considered 
in Fig. [17] 

In addition, the shape of the transition from large to 
small is by no means trivial. One of the key features 
in the radiated energy is that, especially in the P = lP qc 
sequence but also for P = 2P qc , there appear additional 
local extrema which match the number of encounters. 
The observed structure in E rac i shows a remarkably clean 
periodicity as a function of and should be compared 
with corresponding features in the final spin and mass 
in [23]. We find that these features are determined en- 
tirely by the dynamics during the last encounter. Zoom 
whirl effects in the P = lP qc sequence minimize radi- 
ated energy. We find that the energy is less than that of 
a quasi-circular binary in direct contrast to [91j . We will 
interpret these observation in the next section. 

Looking at our findings presented in Fig.|15|and|17| one 
may wonder why the additional peaks, i.e. additional en- 
counters, are present in the lowest momentum sequence, 
but not in the higher momentum ones. The answer lies 
in the initial binding energy. For the large P cases only 
those evolutions which radiate a lot of energy during the 
first encounter will be bound orbits (dynamical captures) . 
As it turns out, the radius of capture for those evolutions 
is inside the ISCO of a single Schwarzschild black hole 
of the same total mass. The capturing encounter generi- 
cally is a whirl and thus an unstable orbit. After the BHs 
are dynamically captured they will have lost additional 
angular momentum and energy. Thus the pericenter dis- 
tance on their next encounter will be even smaller and 
therefore always end in a merger. 

Another observation is that in the P — 2P qc sequence 
there is a second peak next to the global maximum in 
contrast to the higher P sequences. The peak arises from 
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FIG. 14. Radiated energy as a function of the shooting angle O for the P — lP qc runs. The small insets illustrate the 
corresponding orbital dynamics. One can see that the global maximum does not correspond to a zoom-whirl orbit. Strongest 
zoom- whirl behavior is rather associated with the local minimum in E ra d near O = 48.5°. Compare with Fig. |12| 



contributions during the second (and last) encounter. 
While we do observe similar orbital dynamics also for 
the higher P cases, we however do not see a correspond- 
ing peak. The reason is obvious once one compares the 
GW amplitudes during the capturing first encounter with 
the amplitude during merger, see Fig. [8] and [9] For the 
large P evolutions the mergers on the second encounter 
only have a negligible contribution to the radiated en- 
ergy, but the whirly, capturing encounter dominates the 
energy loss. 



Results for unequal mass runs are shown in Fig. |18| 
According to our findings the scaling of radiated energy 
with mass ratio is eccentricity-dependent. Comparing 
the maxima in E rac i between equal mass and unequal 
mass runs we find that a mass ratio 1:2 still gives a max- 
imal efficiency which is not too far away from the corre- 
sponding equal mass run with the same P/Madm- This 
result is in contradiction to our expectation from quasi- 
circular binaries where E rac i decreases steeply with mass 
ratio. Also our results for mass ratios 1:3 show a simi- 
lar trend suggesting that a such mass ratios still are (in 
the above sense) far away from the geodesic limit. Our 
results suggest further parameter studies to analyze the 
scaling in the eccentric regime along the mass ratio axis. 
Clearly, E rac i is much more sensitive to P rather than v. 



C. New diagnostics 

Many interesting questions about BHB cannot be tack- 
led by just looking at gauge-invariant quantities. In this 
section we suggest new diagnostics that are helpful to 
interpret these spacetimes. 

A first example is the observation in |22j that maxima 
in E ra d coincide with a particular orbital configuration 
at the time of merger. Whenever the angle between the 
tangent vector of the puncture and the separation vector 
D at the time of merger is largest, the radiated energy 
is maximized. Here, we confirm this behavior also for 
P = 2P qc orbits in Fig. 19 demonstrating the robustness 



of our gauge-dependent conclusions in [22] . 

We interpret this empirical finding in the following way. 
Maximizing the above mentioned angle translates into 
maximizing L = D x P, the Newtonian expression for the 
angular momentum of two point masses. We therefore 
conjecture based on our data set that the strongest ring- 
down signals are caused by those evolutions which max- 
imize the angular momentum at the moment of merger. 

Another useful diagnostic is the histogram of D(t) (see 
Fig. 20 and [22]). It measures the time the binary spends 
within an interval D ± AD of coordinate separation. We 
focus on two important conclusions drawn from this plot. 
First, as already reported in |22j the whirls show up as 
a sharp and well-defined peak allowing us to measure 
the radii of unstable circular orbits in these highly non- 
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P — lPq C case there are multiple extrema we find only two 
maxima here. Note the fine sampling around the extrema. 
The rather large gap at B = Qbu reflects the problem visi- 
ble in the upper right inset and in Fig. |22| namely that the 
BHs zoom out to very large distances, which implies large t m 
and consequently high computational costs (also due to the 
requirement of higher resolution). 



linear spacetimes. Second, the whirl radii are becom- 



ing systematically tighter as P increases (see Fig. 20) 



which is related to a higher Kerr-parameter of the merger 
remnant HM]. We checked the coordinate separation 
as a function of time separately to exclude a possible issue 
with our merger time estimate which is not well suited 
for large P. 

For low momenta [5TJ[23] showed that the final spin pa- 
rameter lies within 0.6 < a < 0.0832 with the tendency 
that the spin parameter grows with the initial momen- 
tum. Thus the resulting background spacetime will have 
a tighter ISCO. As the binary spends considerably more 
time at the whirl radius than at a Newtonian pericenter 



at the same distance (see Fig. 20) the binary radiates 
much more efficiently if this whirl occurs at a smaller ra- 
dius. Thus the geodesic analog together with our gauge- 
dependent diagnostics give a natural explanation for our 
earlier observations: In the high momentum case zoom- 
whirl orbits do coincide with the most efficient radiators 
where the whirls are tight, while this is not the case in 
the low momentum regime, where the whirl radii are sig- 
nificantly larger. 

As a final diagnostic we present trajectories of the bina- 



ries through phase space, see Figs.[2TJ[22j[23j and|24| We 

choose D{t) and the coordinate velocity v(t) of the punc- 
tures as generalized coordinates. This construction is ex- 
plicitly coordinate dependent, and switching to another 
gauge will lead to different trajectories. However, previ- 
ous investigations led to the conclusion that the moving 
puncture gauge leads to puncture tracks that correspond 
rather well to what an observer sees from infinity, e.g. 



FIG. 16. Radiated energy for P = 10P qc (i.e. P « 0.6M). 
The values for O are not comparable with those from the 
other P-sequences as the initial separation was chosen to be 
D(t = 0) = 50M due to the large junk radiation. Only at 
P = WPqc are we able to exceed the Hawking limit [101] on 
the energy release of two Schwarzschild BHs far apart without 
orbital angular momentum. With P = 10P qc the spacetime 
is very different from Schwarzschild. Note that the results are 
below the current maximum reported value of 35 ± 5%. The 
contribution from junk radiation is about 0.05 and is included 
in the data shown. 



[?2l 15B1 1103] . In particular, for orbiting motion one can 
argue based on the shift condition that this should be the 
case |SS], although for linear motion the situation is dif- 
ferent. Hence, one has to keep the gauge issue in mind, 
but the moving puncture gauge leads to rather robust 
features in the phase space trajectories as we will discuss 
next. 

To familiarize oneself with the trajectories in phase 
space consider a circular motion with constant velocity. 
This motion corresponds to a single point in a D-v phase 
space diagram. A Kepler ellipse corresponds to a line, 
which is curved according to Kepler's third law. 

In Fig. [21] we show runs for P = lP qc for four different 



0, while Fig. 22 gives a global impression of many differ- 
ent angles for P = 2P qc (each for equal masses). These 
simulations start at the lower right corner at D = 20M, 
v = and quickly rise due to the initial gauge adjust- 
ment, so that the coordinate velocity of each BH ap- 
proaches the value given expected from the intial data. 
Then the trajectory moves towards the upper left as the 
orbits shrink and the punctures move faster. When the 
black holes merge the trajectory ends in the origin at 
D = 0, v = 0. Whirls or parts of tight circular orbits are 
indicated by approximately constant D but decreasing 
v, with D s» 2M-5M. Zooms follow roughly the shape 
of elliptic orbits, with D varying between 5M-10M at 
pericenter out to an apocenter at 15M-27M in Fig. [2l] 
or up to 40M in Fig. [22j Fig. [22] also shows one orbit for 
P = 2P qc which escapes to infinity. 
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FIG. 17. Radiated energy for the lqc-6qc runs (i.e. 0.06M > P > 0.36A4") including the values corresponding to a quasi- 
circular binary and the Hawking limit. The vertical lines mark Qbu for each momentum scale considered. This value is close 
but not identical to the threshold of immediate merger. As expected the shooting angle where zoom-whirls occur is closer to 
the = (head-on) case for higher initial momentum. The higher the initial momentum the more energy is radiated. In the 
lP qc sequence we have not included the data for Q = 60°, but the efficiency of radiation agrees within plotting accuracy with 
the value for a quasi-circular BHB as expected from extrapolating the data set shown. 



A head-on collision in phase space looks very much 
like an eccentric binary starting on the v = line at the 
value of D that corresponds to the same total energy. We 
note that a head-on collision always constitutes an upper 
envelope in the phase space, i.e. fixing P and comparing 
= runs with 0^0 evolutions we always find vit, = 
0) > v(t, 9^0) (at least for D(t) < D(t = 0)). 

We note that the motion of the punctures at the onset 
of merger is still rather mildly relativistic. This find- 
ing turns out to be surprisingly insensitive to the initial 
momentum. Fig. [23] shows runs for different P with an- 
gles chosen for maximum radiation efficency. Increasing 
the initial momentum of the punctures leads to a mo- 
tion which is rather relativistic when entering the whirl 
phase, up to v ~ 0.5 and Lorentz factor W ~ 1.155 for 
P = 6P qc (see Fig. 



23 ) , but in the whirl they decelerate 



by large amounts (for" P = 6P qc the decrease in velocity 
is as large as during the merger). At merger however 
all equal mass evolutions approach v = 0.22 ± 0.02 or 
W = 1.025 ± 0.005. 

For a detailed look at the transition zone between 



bound and unbound evolutions, see Fig. 22 for the P — 
2P qc sequence varying 0. First, we recognize the two ex- 
treme cases of merger and unbound motion. Inbetween 
we find a very complex transition which is governed by 



several additional (non-closed) loops corresponding to or- 
bits that zoom out after a first whirl phase thereby slow- 
ing down (i.e. move to the lower right) before returning 
towards the upper left. Note, that during the next ap- 
proach the binary must follow a path further to the lower 
left because the system is dissipative. The set of apocen- 
ters from all runs in this sequence forms a lower envelope 
(as the upper envelope mentioned before) that is never 
crossed by any of our evolutions at the same initial D. 
Note how far the zooms may extend when the binary 
approaches the merger / fly-by threshold resulting in an 
ever larger runtime. We face difficulties evolving orbits 



near — s- 



merger 
fly-by 



because those orbits need very long 



evolutions (decrease of accuracy) and the black holes may 
reach distances close to or even beyond the wave extrac- 
tion sphere, thereby producing artificial features in E ra d- 

An investigation of the phase space trajectories of un- 
equal mass BHBs reveals that the general shape of the 
trajectories is quite robust with respect to the mass ra- 
tio. An example is given in Fig. [24] There are two basic 
differences: (1) Unequal mass orbits move systematically 
slower compared to a corresponding equal mass binary 
with comparable initial momentum. (2) Unequal mass 
mergers start their final plunge from an inreasingly larger 
coordinate separation than equal mass mergers. With re- 
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Efficiency of gravitational radiation 




Shooting angle e[°] Coordinate Separation D[M] 



FIG. 18. Radiated energy E ra d for mass ratios 1:1, 1:2 and 
1:3 and some particular values of P. Qualitatively the re- 
sults for unequal masses are similar to those for equal masses. 
For example, note the drop in E ra d after a global maximum, 
the global maximum corresponds to an orbit that roughly 
completes 1 orbit, and zoom- whirl behavior for low P is asso- 
ciated with inefficient radiation as in the equal mass regime. 
For P = 2P qc , zoom-whirl behavior for mass ratio 1:2 occurs 
for larger Q than in the 1:1 case. 



FIG. 20. Histogram of the coordinate separation D for three 
different evolutions with P — 2,A,QP qc . This plot shows how 
much time a binary has spent at a given separation D. All 
runs shown here correspond to the longest whirl phase found 
at each fixed momentum. Clearly visible is the tightening of 
the whirl radius for larger P (compare with [22]). The overlap 
of the shaded region in P = 6P qc with the histogram happens 
because during the whirl the separation is indeed shorter than 
at the onset of merger. 




-4 -2 2 4 



X[M] 

FIG. 19. The inner region of puncture tracks from two 
different runs. The P = 2P qc , O = 24.1° run (blue solid 
line) corresponds to the global maximum in E ra d. The yel- 
low straight lines through the origin represent the separation 
vectors D at the time when a common horizon forms. The 
result [33J that most efficient mergers occur when the tangent 
vectors of the orbits are closest to being orthogonal to D(t m ) 
carries over to larger P. 
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FIG. 21. Phase space for P = lP qc sequence. The 9 = 48.2 
evolution takes a detour in phase space thereby avoiding the 
region where radiation is most efficient. The O = 47 evolution 
radiates more efficiently because it reaches further towards the 
upper left and at the same time spends considerable time at 
low D (see Fig. 6 in [22]). 



gard to gravitational radiation, unequal mass binaries do 
not tap as deeply into the gravitational potential as equal 
mass binaries, thus cannot extract as much energy from 
the spacetime. Furthermore, the deceleration is milder 



for higher mass ratios, again suggesting weaker signals 
from higher mass ratios. 
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FIG. 22. Phase space for P — 2P qc sequence. All runs except 
B = 25.5° are bound. Near the threshold of immediate merger 
the zooms (24.3° < O < 24.9°) may extend far out. The inset 
shows a close-up of the inner region where whirls occur (the 
point density is enlarged there). 



Phase space (momentum dependence) 





-- ip, c ,e=47° 




.... 2P 5C ,9=23.9° 




4P 5c ,e=15.3° 




• ■ 6P,„e=13.4° 


ei*7"-». 

f: 




_f * - * 




I 

; 





0.0 1 ' 1 ' 1 

5 10 15 20 

Coordinate separation 



FIG. 23. Phase space for different initial momenta P with 
angles chosen for near maximal E ra d- Interestingly, close to 
the merger all binaries reach the same coordinate velocity 
independent of the large differences in initial conditions. This 
effect appears as a blurred "focal" point in phase space. Hence 
the whirl itself becomes more important for E ra d than the 
actual merger. The deceleration in the 6P qc whirl phase is 
larger than during the merger. 



IV. CONCLUSIONS 

Numerical relativity has confirmed the existence of 
zoom-whirl orbits beyond the geodesic and PN regime 
[EH221 EH [26] , thereby emphasizing their general rela- 
tivistic origin. Previous studies explored the rich exten- 
sion to the phenomenology of the GR two-body problem 



FIG. 24. Mass ratio 1:3. The whirl phase is slightly but 
not significantly longer than for the equal mass configuration. 
This difference can be understood by comparing radiation ef- 
ficiencies of these systems. A comparison of the 22-mode to 
the equal mass runs suggests strong sensitivity on the mass 
ratio. Note that the merger is not the strongest signal in ^4. 



offered by zoom-whirl dynamics in various ways. In this 
work we performed numerical relativity simulations to in- 
vestigate the parameter space of comparable mass, non- 
spinning, eccentric BHB for low and intermediate mo- 
menta more comprehensively than before. We explored 
zoom-whirl behavior in both the hyperbolic and elliptic 
regime carrying out more than 100 numerical evolutions 
in order to obtain a decent sampling of the underlying pa- 
rameter space. We discussed various features of the orbits 
and characterized the corresponding GW emission, and 
we developed new diagnostics to analyze binary space- 
times by using phase space trajectories and a histogram 
of the coordinate separation. 

For elliptic orbits, we discover zoom-whirls with im- 
prints in the GWs that are comparable in amplitude to 
the merger waveform for eccentricities as low as e <~ 0.5. 
This is an important finding for the astrophysical rele- 
vance of zoom-whirl orbits. In particular, such values 
are within expected eccentricities of supermassive BHB 
that have resulted from galaxy mergers and subsequent 
star- or gas-driven hardening |49j . For low momenta, 
there occur several minima and maxima in the radiated 
energy when varying the shooting angle from head-on to 
quasi-circular orbits. We demonstrate that zoom-whirl 
dynamics may actually minimize the radiated energy in 
sharp contrast to [32] . 

In the elliptic regime whirls are found only during their 
last encounter. They emerge in disjunct intervals of the 
initial angular momentum (i.e. shooting angle). Appar- 
ently, as long as the binaries are not circularized just prior 
to plunge, zoom-whirls can always be found during the 
last encounter by a very modest amount of fine-tuning. 
In the hyperbolic regime we find that all evolutions that 
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lead to dynamical capture reveal whirl features during 
the capturing encounter and then simply plunge during 
the following encounter, potentially from a large separa- 
tion. 

High-momentum zoom-whirls maximize the radiated 
energy. The Kerr spacetime that the merger remnant 
will settle down to exhibits a larger spin parameter. This 
translates into a tighter unstable circular orbit (i.e. the 
whirl radius) resulting in more pronounced dynamics of 
the mass quadrupole. Especially the first, capturing 
close-encounter burst in high momentum evolutions can 
easily overwhelm the merger signal. 

In the unequal mass case eccentric BHB are found 
to be more efficient radiators than expected from quasi- 
cirular studies of unequal mass BHB. As a consequence 
the whirls are not significantly longer for the mass ratios 
under consideration. We note that numerical relativity 
is seeing improvements in dealing with large mass ratios 
[1041 1105] . More detailed studies for unequal masses and 
also the inclusion of spin are promising directions for ec- 
centric BHB simulations in the future. 



The present work as well as other studies strongly sug- 
gest to include eccentricity in the waveform templates 
used in the data analysis of GW detectors, since eccen- 
tricity effectively breaks degeneracies in parameter space. 



ACKNOWLEDGMENTS 

We are grateful to David Hilditch, Sebastiano 
Bcrnuzzi, Ulrich Sperhake, Frans Pretorius, Luis Lehner, 
Deirdre Shoemaker, Doreen Miiller, Norbert Lages, and 
Marcus Thicrfeldcr for discussions, and we also thank 
Marcus Thierfeldcr for giving access to his event hori- 
zon finder. This work was supported in part by grants 
DFG GK 1523, DFG SFB/Transregio 7, and DLR LISA 
Germany at the Friedrich-Schiller University Jena as 
well as'NSF Grant PHY-0963136, and NASA Grant 
NNX11AE11G at the University of Illinois at Urbana- 
Champaign. Computations were performed at LRZ Mu- 
nich. 



[3] 
[4] 



[5] 

[6] 
[7] 



K. Glampedakis and D. Kennefick, Phys. Rev. D66, 

044002 (2002), |arXiv:gr-qc/0203086| 

The term zoom-whirl is attributed to C. Cutler and E. 

Poisson (who in turn mention it was suggested by K. 

Thorne). 

R. A. Hulse and J. H. Taylor, Ap,] 195, L51 (1975) 
M. Burgay, N. DAmico, A. Possenti, R. Manch- 
ester, A. Lyne, B. Joshi, M. McLaughlin, M. Kramer, 
J. Sarkissian, F. Camilo, V. Kalogera, C. Kim, and 
D. Lorimer, [Nature| 426, 531 (2003), |astro-ph/031207l"1 
M. Kramer, I. Stairs, R. Manchester, M. McLaugh- 
lin, A. Lyne, R. Ferdman, M. Burgay, D. Lorimer, 
A. Possenti, N. DAmico, J. Sarkissian, G. Hobbs, 
J. Reynolds, P. Freire, and F. Camilo, Science 314, 97 
(2006), |arXiv:astro-ph/0609417 [astro-ph]| 
A. Sillanpaa, S. Haarala, M. Valtonen, B. Sundelius, 
and G. Byrd, Astrophys.J. 325, 628 (1988) 
M. Valtonen, S. Mikkola, D. Mer ritt, A. Gopaku- 
mar, H. Lehto, et al, |Astrophys.J.| 709, 725 (2010), 
|arXiv:0912.1209 [astro-ph.HEj | 
J. 



28, 



[101 

[11] 
[12] 
[13] 
[14] 
[15] 



Rev. D77, 103005 



Levin and H. Contreras, Class. Quant. Grav 
175001 (2011), |arXTv:1009 2533 [gr-qcj| 
S. Chandrasekhar, The Mathematical Theory of Black 
Holes (Oxford University Press, USA, 1983) 
K. Martel, Phys. R ev. D69, 044025 (2004), arXiv:gr 
qc/0311017 [gr-qc]| 
J. Levin and G. Perez-Giz, Phys 

(2008) , |arXiv:0802 0459 [gr-qc]| 

S. Dra sco and S. A. HughesTPhys. Rev. D73, 024027 
(2006), |arXiv:gr-qc/0509101 [gr-qc]| 
R. Haas, |Phys. Rev.| D75, 
|arXiv:0704.0797 [gr-qcjp 
J. Levin and B. Grossman, Phys. 

(2009) , [arXiv:0809.3838 [gr-qc]| 
R. Grossman and J. Levin, Phys. 
(2009),|arXTv:0811.3798 [gr-qc] 



[16 

[I?: 
[is: 

[19 
[20 
[21 
[22 

[23; 



124011 (2007), 



Rev. D79, 043016 



Rev. D79, 043017 



24 

[25 
[26; 

[27 

[28; 



[29 
[30 



S. A. Hughes, S. Drasco, E. E. Flanagan, and 
J. Franklin, |Phys. Rev. Lettl. 94, 221101 (Jun 2005) 
M. E. Pati and C. M. Will, Phys. Rev. D65, 104008 
(2002) 

F. Pretorius and D. Khurana, Cla ss. Quant. Grav.| 24, 
S83 (2007), |"arXiv:gr-qc/0702084| 

U. Sperhake, E. Berti, V. Cardoso, J. A. Gonzalez, 
B. Briigm ann, and M. Ansorg, |Phys. Rev.| D78, 064069 
(2008), | arXiv:0710.3823 [gr-qc]| " 

I. Hinder, B. Vaishnav, F. Herrmann, D. Shoemaker, 
and P. Laguna, |Phys. Rev.| D77, 081502 (2008), 
|arXiv:0710 5167 [gr-qcj| " 

M. C. Washik, J. Healy, F. Herrmann, I. Hinder, D. M. 
Shoemaker, P. Laguna, and R. A. Matzner, Phys. Rev. 
Lett. 101, 061102 (2008), |arXiv:0802.2520 [gr-qc]| 
R. Gold and B. Briigmann, Class. Quant. Grav. 27, 
084035 (2010), [arXTv:0911 3862 [gr-qc]| 
J. Healy, P. Laguna, R. A. Matzner, and D. M. Shoe- 
make^ Phys. Rev. D81, 081501 (2010), [arXiv:0905.3914| 
[gr-qc] 

J. Healy, J. Levin, a nd D. Shoemaker, |Phys. Rev. Lett.| 
103, 131101 (2009),[arXTv:0907.0671 [gr-qcj| 



C. Cutler, D. Kennefick, and E. Poisson, |Phys. Rev. D 
50, 3816 (Sep 1994) 

J. Healy, F. Herrmann, I. Hinder, D. M. Shoemaker, 
P. Laguna, and R. A. Matzner, Phys. Rev. Lett. 102, 
041101 (2009), [arXiv:0807.3292 [gr-qc]| 
K. Martel and E. Poisson, Phys. Rev. D60, 124008 
(1999) 



K. Martel, AIP Conf. Proc. 493, 48 (1999), arXiv:gr- 
qc/9908043 [gr-qc] | 

D. A. Brown and P. J. Zimmerman, Phys. Rev. D81 
024007 (2010), |arXiv:0909.0066 [gr-qc]] 
T. Cokelaer and D. Pathak, Class. Quant. Grav. 26. 
045013 (2009), [arXTv:0903 .4791 [gr-qc]| 



18 



B. Vaishnav, I. Hinder, D. Shoemaker, and F. Herrman, 
Classical and Quantum Gravity 26, 204008 (2009) 
B. Mikoczi, B. Kocsis, P . Forgacs, and M. Vasuth(2012), 
|arXiv:1206.5786 [gr-qc]| 



[astro-ph.HE] 

W. E . East, F. Pretorius, an d B. C. S tephens , |Phys. 
RevT| D85, 124009 (2012), |arXiv:1111.3055 [ astro- 
ph.HE] 

R. Gold, S. Bernuzzi, M. Thierfelder, B. Briigmann, and 
F. Pretorius (20 11), |arXiv: 1109.5128 [gr^qcj] 
B. C. Stephens, W. E. East, and F. Pretorius, Astro- 
phys. J. 737, L5 (2011), |arX"iv:1105.3175 [astro-ph.HEj 
M. Shibata, H. Oka wa, and T. Ya mamoto, Phys. Rev. 



D78, 101501 (2008), |ajXTv:0810.47 35 [gr-qc] 
U. Sperhake, V. Cardoso, F. Pretorius, E. Berti, T. Hin- 
derer, et al, |Phys. Rev. Lett.| 103, 131102 (2009), 
|arXiv:0907.1252 [gr-qc] p 

U. Sperhake, V. Cardoso, F. Pretorius, E. Berti, and 



J. A. Gonzalez. 
|arXiv:0806.l738 

E. Berti, V. Cardoso 

F. Pretorius, et al. 



Phys. Rev. Lett, 
gr-qc] 



101, 161101 (2008), 



Rev 



Hinderer 
D81 



M. 
104048 



Lemos, 
(2010), 



598, 419 (2003), arXiv:astro- 



Anderson, H. Baumgardt, 
591, L131 (2003), 

Hamilton, Astro- 



[59] 



L. Wen 

ph/0211492 [astro-phj| ' 
D. Pooley, W. Lewin 
A. Filippenko, et al, Astrophys.J 
arXiv:astro-ph/0305003 [astro-ph] 
M. C. Miller and D. P. 
phys. J. Lett. (2002), |arXiv:astro-ph/0202298 [astro-ph] 
W. H. Lee, E. Ramirez-Ruiz, and G. van de Ven, Astro 
phys. J.| 720, 953 (2010), [arXiv:0909. 2884 [astro-p h.HE] 
R. M. O'Leary, B. Ko csis, and A. Loeb(2008), 

|arXiv:0807.2638 [astro-ph] | 

W. R. Brown(2008), |arXiv:0811.0571 [astro-ph]] 
S. F. Portegies Zwart, H. Baumgardt, S. L. McMil- 
lan, J. Makino, P. Hut, et al, Astrophys.J. (2005), 
arXiv:astro-ph/0511397 [ast ro-ph]| 

O. Blaes, M. H. Lee, and A. Socrates, Astrophys.J. 578, 
775 (2002), |arXiv:astro-ph/0203370 [astro-ph]| 



[60 
[61 



W. E. East a nd F. Pretorius(2012), |arXiv:1208.5279| 



Phys 

|arXiv:1003.0812 [gr-qcj| " 

U. Sperhake, V. Cardoso, F. Pretorius, E. Berti, T. Hin- 
derer, et al.(2010), [arXiv: 1003.0882 [gr-qc]| 
P. C. Peters, Phys. Rev. 136, B1224 (1964) 
A. Sesana, C. Roedi g, M. Reyno lds, and M. Dotti(2011), 
|arXiv:1107.2927 [astro-ph.CO]| 

F. Kha n, A. Just, and D. Merritt,|Astrophys.J.] 732, 89 

(2011) , [arXiv:1103.0272 [astro-ph. CO]T 
M. Preto, I. Berentzen, P. Berczik, and 
R. Spurzem(2011), |arXiv:1102.4855 [astro-ph.GA]| 

J. C. B. Papaloizou, R. P. Nelson, and F. Masset, |As-| 

tronomy and Astrophysics 366, 263 (jan 2001) 

P. J. Armitage and P. Natarajan, Astrophys.J. 634, 921 

(2005), |arXiv:astro-ph/0508493 [astro-phj| 

J. Cuadra, P. Armitage, R. Alexander, and M. Begel- 

man(2008), |arXiv:0809.0311 [astrc^ph]1 

C. Roe dig and A. Sesana, |J.Phys.Conf.SerTl 363, 012035 

(2012) , [arXivillll. 3742 [astro-ph.C^jp 

S. F. Portegies Zwart and ~ S~~McMillan(1999), 
arXiv: astro-ph /9910061 [astro-ph] 

M. Benacquista, S. F. Portegies Zwart, and F. Ra- 
sio, |Class.Quant.GravT| 18, 4025 (2001), |arXiv:gr-| 
qc/0010020 [gr-qc] 

Astrophys.J. 



[62 
[63 
[64 
[65 

[66 
[67 

[6s; 

[69 

[70; 

[71 

[72; 



[73 



"[74" 

[75; 
[7e; 
[77; 
[7s; 



[79; 

[80 
[81 

[82 
[83; 



[85 



C. Hopman and T. Alexander, Astrophys.J. 629, 362 
(2005), arXiv:astro-ph/0503672 [astro-ph] 
I. Mandel, D. A. Brown, J. R. Gair, and M. C. Miller, 
Astrophys.J. 681, 1431 (2008), arXiv:0705.0285 
Ph] 



astro- 



M. Iwasawa, Y. Funato, and J. Makino, [As trophys.J. 
651, 1059 (2006), |arXiv:astro-ph/0511391 |astro-ph|| 
M. Volonteri, F. Haardt, and P. Madau, |Astrophys.J.| 
582, 559 (2003) , |arXiv:astro -ph/02072 76 [astro-ph] | 
F. Antonini and H. Perets,|Astrophys.J.| 757727~(2012) , 
|arXiv:1203.2938 [astro-ph. GAj| 

I. Kowalska, T. Bulik, K. Belczynski, M. Dominik, 
and D. Gondek-Rosinska, Astron.Astrophys. 527, A70 
(2011), |arXiv:1010. 0511 [astro-ph. CO]] 
B. Kocsis, A. Ray, and ST P. Zwart(2011), 
|arXiv:1110.6172 [astro-ph. G A] | 

B. Briigmann, J. A. Gonzalez, M. Hannam, S. Husa, 
U. Sperha ke, and W. Tichy |Phys. Rev.| D77, 024027 
(2008), |arXiv:gr-qc/ 0610128 

B. Briigmann, W. Tichy, and N. Jansen, Phys. Rev. 
Lett. 92, 211101 (2004), [gFqc/0312112| 
S. Husa, J. A. Gonzalez, M. Hannam, B. Briigmann, and 
U. Sperhake, Class. Qu ant. Grav. 25, 105006 (2008), 
|arXiv:0706.0740 [gr-qc]| 

S. Brandt and B. Briigmann, Phys. Rev. Lett. 78, 3606 
(1997), [gr^qc/9703066| 

M. Ansorg, B. Briigmann, and W. Tichy, Phys. Rev. 

D70, 064011 (2004), [gpqc/0404056| 

M. Campanelli, C. O. Lousto, P. Marronetti, and 

Y. Zlochower, Phys. Rev. Lett. 96, 111101 (2006), gr- 

qc/0511048| 

J. G. Baker, J. Centrella, D.-I. Choi, M. Koppitz, and 
J. van Meter, Phys. Rev. Lett. 96, 111102 (2006), gr- 
qc/0511103| 



M. Shibata and T. Nakamura, Phys. Rev. D52, 5428 
(1995) 

T. W. Baumgarte and S. L. Shapiro, Phys. Rev. D59, 
024007 (1998), gr-qc/9810065 

E. Schnetter, IClass. Quant. Grav.| 27, 167001 (2010), 
|arXiv: 1003.0859 [gr-qcjp 

B. Walther, B. Briigmann, and D. Miiller, Phys. Rev. 
D79, 124040 (2009), [arXTv:0901.0993 [gr-qc]| 
Generally speaking, if such an orbit is unbound then the 
black holes separate without close encounter. If bound, 
then after reaching the first apocenter the orbit becomes 
inward-bound, which in principle is already included for 
6 6 [0,90°]. 

W. H. Press, B. P. Flannery, S. A. Teukolsky, and W. T. 
Vetterling, Numerical Recipes in C, 2nd ed. (Cambridge 
University Press, New York, 1992) 

R.-M. Memmesheimer and G. Schafer, Phys. Rev. D71, 



044021 (2005), arXiv:gr-qc/0411057 [gr-qc 
D. Miiller and B. Briigmann, Class. Quant. Grav. 27, 
114008 (2010), |arXiv:0912. 3125 [gr-qcj| 
D. Miiller, J. Grigsby, and B. Briigmann, Phys. Rev. 
D82, 064004 (2010), [arXiv: 1003.4681 [gr-qc]| 
D. Alic, L. Rezzolla, I. Hinder, and P. Mosta, 



Class.Quan t.Grav.| 27, 245023 (2010), |arXiv:1008.2212| 
[gr-qc] 

M. Thierfelder, "Event horizon finder for black holes," 

(2008), Diploma Thesis, University of Jena 

M. Thierfelder and B. Briigmann, Phys. Rev. D (to be 

submitted) 



19 



M. Hannam, S. Husa, D. Pollney, B. Briigmann, and 
N. 6 Murchadha, Phys. Rev. Lett. 99, 241102 (2007), 
|gr-qc/0606099| 

B. Rei mann and B. Bri igmann, Phys. Rev. D69, 124009 
(2004), |gFqc/0401098| 

M. Alcubierre, Introduction to 3+1 Numerical Relativity 
(Oxford University Press, USA, 2008) 
E. Berti, V. Cardoso, J. A. Gonzalez, U. Sperhake, and 
B. Briigmann, |Class. Quant. GravT] 25, 114035 (2008), 
|arXiv:071 1 . 1097 [gr-qcjp 

S. Shapiro and S. Teukolsky, Blackholes, White Dwarfs 
and Neutron Stars: The physics of compact objects 
(John Wiley & Sons, 1983) 

P. C. Peters and J. Mathews, P hys. Rev. 131, 4 35 (1963) 
M. Tes smer and G. Schafer, |Annalen Phys. 523, 813 
(2011), |arXiv:1012.3894 [gr-qc]| ~ 

I. Hinder, F. Herrmann, P. Laguna, and D. Shoemaker, 
|Phys. Rev.| D82, 024033 (2010), |arXiv:0806. 1037 [gr-qc]| 
J. Hough (1994), prepared for Edoardo Amaldi Meeting 
on Gravitational Wave Experiments, Rome, Italy, 14-17 
Jun 1994 

S. Babak, M. Hannam, S. Husa, and B. Schutz, "Re- 
solving Super Massive Black Holes with LISA," (2008), 



[97] 



arXiv:0806.1591 [gr-qc] 




P. Amaro-Seoane, S. Aoudia, S. Babak, P. Binetruy, 
E. Berti, et oZ.(2012), arXiv:1201.3621 [astro-ph.CO] 


K. Yagi, Class. Quant. Grav. 


29^ 075005 (2012), 


|arXiv:1202.3512 [astro-ph.CO] 
M. Hannam and I. Hawke, 


Gen.Rel.GravT] 43, 465 


(2011), arXiv:0908.3139 [gr-qc] 





Proa] 1682(2010), 



[100] 

[101] 
[102] 
[103] 



[104] 
[105] 



C. Van Den Broeck, 
|arXiv:1003.1386 [gr-qc] | 
J. S. K ey and N. J. Cornish, |Phys. Rev.| D83, 083001 
(2011), {arXiv:1006.3759 [gr-qc] | 
S. Hawking, Phys. Rev. Lett. 26 , 1344 (1971) 
U. Sperhake, [Lect. Notes Phys.| 769, 125 (2009) 
P. Diener, F. Herrmann, D. Pollney, E. Schnetter, 
E. Seidel, R. Takahashi, J. Thornburg, and J. Ventrella, 
Phys. Rev. Lett. 96, 121101 (200 6), |gr-qc/0512108| 
C. O. Lousto and Y. Zlochower, Phys. Rev. Lett. 106, 
041101 (2011), [arXiv:1009. 0292 [gr-qc] | 
J. A. Gonzalez, U. Sperhake, and B. Briigmann, Phys. 
Tlev]D79, 124006 (2009), [arXTv:0811.3952 [gr-qc 



